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The majority of magnetoelectric crystals show complex temperature-magnetic field or 
temperature-pressure phase diagrams with alternating antiferromagnetic incommensurate, magne- 
toelectric, and commensurate phases. Such phase diagrams occur as a result of successive magnetic 
instabilities with respect to different order parameters, which usually transform according to differ- 
ent irreducible representations (IR) of the space group of the crystal. Therefore, in order to build 
a phenomenological theory of phase transitions in such magnetoelectrics one has to employ several 
order parameters and assume the proximity of various instabilities on the thermodynamic path. 
In this work we analyze the magnetoelectrics MnWC>4, CuO, NaFeSi2C>6, NaFeGeaOe, Cu3Nb20s, 
Q-CaCr204, and FeTe20sBr using the praphase concept and the symmetry of the exchange Hamil- 
tonian. We find that in all the considered cases the appearing magnetic structures are described by 
IR's entering into a single exchange multiplet, whereas in the cases of MnWCU and CuO by a single 
IR of the space group of the praphase structure. Therefore, one can interpret the complex phase 
diagrams of magnetoelectrics as induced by a single IR either of the praphase or of the symmetry 
group of the exchange Hamiltonian. Detailed temperature-magnetic field phase diagrams of MnWCU 
and CuO for certain field directions are obtained and the magnetic structures of the field-induced 
phases are determined. 

PACS numbers: 75.85.+t, 77.84.-s 



I. INTRODUCTION 



In the recent decade magnetoelectricity has become 
one of the focal points of interest for both the mag- 
netic and ferroelectric communities. The potentialities 
of combining magnetic and ferroelectric properties in 
the same system and, which is more important, the 
possibility of controlling one property by altering the 
other, open the way to new interesting applications. The 
promising applications of the magnetoelectric effect in- 
clude new types of devices for reading, writing, and stor- 
age of information^^ various sensorsj&i microwave^ and 
spintronics 6 ' devices, and wireless energy transfer and en- 
ergy harvesting technologies.- Therefore considerable ef- 
forts are devoted from both the theoretical and experi- 
mental points of view to the search of new magnetoelec- 
tric crystals and explanation of their properties. 

The vast experimental data accumulated during the 
past years allows establishing specific features pertinent 
to magnetoelectrics, which have been summarized in re- 
cent reviews. The magnetoelectrics can be divided into 
two classes^ To the so-called type-I magnetoelectrics, 
which are sometimes referred to as ferroelectromagnets, 
belong crystals in which ferroelectricity and magnetic or- 
der appear independently and have different sources. The 
prominent examples are BiFe03- and YMn03^ Type-I 
multifcrroics are usually characterized by rather high fer- 
roelectric transition temperatures and large electric po- 
larization (~10 - 100 ^C/cm 2 ), but the generally large 
difference between the ferroelectric and magnetic tran- 
sition temperatures and the different causes of the two 
orders result in small coupling between them^ 



The second class, the type-II magnetoelectrics, is com- 
prised of multiferroics, in which electric polarization ap- 
pears upon magnetic phase transitions. This class is gen- 
erally characterized by much lower transition tempera- 
tures (10 - 40 K) and rather low electric polarization val- 
ues (usually of the order of 10 - 100 /iC/m 2 ). The promi- 
nent examples of such magnetoelectrics are rare-earth 
manganites RMn03 (R=Gd, Tb, and Dy). 11 Compared 
to ferroelectromagnets the type-II multiferroics provide 
much stronger magnetoelectric coupling due to the fact 
that electric polarization is induced by magnetic order. 
Therefore, application of external magnetic field alters 
the magnetic structure, which may result in suppression, 
induction, rotation, or a flop of electric polarization. In 
this paper we focus on the type-II magnetoelectrics. 

Two main features of the type-II magnetoelectrics 
can be established from the analysis of the vast ex- 
perimental data available to date. The first one con- 
sists in the fact that magnetoelectrics often exhibit 
complex temperature-pressure or temperature-magnetic 
field phase diagrams with alternating magnetically or- 
dered modulated, commensurate, and magnetoelectric 
phases. Common appearance of incommensurate mag- 
netically ordered phases and the frequent emergence of 
electric polarization in them determine the description of 
magnetoelectricity from both the micro- and macroscopic 
points of view. Despite the fact that ferroelectricity can 
be also induced by commensurate magnetism (as, for 
example, in rare-earth manganates RM^Os, Ref. 12]), 
this manifests itself in consideration of various complex 
magnetic orders such as screw, cycloidal, helix, and oth- 
ers as sources of magnetoelectricityj 13 ' 14 

From the macroscopic crystal symmetry point of view 



the close connection between the appearance of modu- 
lated and ferroelectric phases in magnetoelectrics was 
recently pointed out by the present authors. 15 It was 
pointed out, in particular, that in some magnetoelectrics 
the transformational properties of magnetic order pa- 
rameters are described by the irreducible representations 
(IR) not satisfying the Lifshitz criterion. This results 
in long-periodical modulation of the magnetic order, as 
well as in ferroelectric states among the low symmetry 
phases induced by these order parameters. Unlike other 
phenomenological models (see, for example, Refs. [la - 
Il8l ] ) , in which one employs separate order parameters for 
different experimentally observed modulation wave vec- 
tors, our approach provides the whole picture description 
using the minimal number of order parameters. In the 
case when a commensurate magnetically ordered phase 
is present in the phase diagram, the order parameters 
belonging to its wave vector are chosen and the incom- 
mensurate phases are found to be well described by the 
existing Lifshitz invariants. This is, for example, the case 
of RMnaOsji 5 - MnW0 4 ^ and CuO 3 In other cases 
(such as, for example, in CuCk), Ref. [21|), when the mag- 
netic order does not lock-in to a commensurate structure, 
one may choose the commensurate wave vector closest to 
that of the incommensurate phases to define the order 
parameters. This approach is justified by the fact that 
such models as, for example, the ANNNI model, which 
describe the magnetic structure by the Ising or Heisen- 
berg Hamiltonians, show a multitude of incommensu- 
rate phases, but nevertheless usually have commensurate 
magnetic structure as the ground state i 22 ' 23 

Another well established peculiarity of magneto- 
electrics consists in the following. Neutron diffraction 
studies of many magnetoelectrics and the respective sym- 
metry analysis reveal that the description of the exper- 
imentally observed magnetic structures requires two or 
more magnetic order parameters belonging to the same 
wave vector . 16 i 24 ~ — These order parameters transform 
according to a single or different IR's of the symmetry 
group of the paraelectric and paramagnetic phase (para- 
phase) and induce electric polarization acting simultane- 
ously. Indeed, in the case where the space group of the 
initial paraphase possesses spatial inversion and the re- 
sulting magnetic cell of the magnetically ordered phase 
coincides with the initial chemical cell (i.e., when k = 0) 
magnetic order parameters transforming according to a 
single IR cannot induce polarization.— However, this is 
generally not the case for fc 7^ 0, for which a single mag- 
netic order parameter may induce electric polarization* 1 ^ 

It is found that, for example, in TbMn03 24 and 
MnWO^ two different order parameters condense suc- 
cessively, whereas in CuCl a 21 ' 28 they condense simulta- 
neously. These order parameters transform according to 
different irreducible representations of the crystal space 
group. Therefore, for the correct phenomenological de- 
scription of the magnetic phase transitions in these mag- 
netoelectrics it is necessary to assume the proximity of 
instabilities with respect to different irreducible represen- 



tations on the thermodynamic path i 18 ' 29 

In this paper we analyze several magnetoelectrics 
(MnW0 4 , CuO, pyroxenes NaFeSi 2 6 and NaFeGe 2 6 , 
Cu3Nb208, a-CaCr204, and FeTe20sBr) from the point 
of view of crystal and exchange symmetry. We show that 
in all the considered cases their magnetic structures are 
described by IR's entering into a single exchange mul- 
tiplet. For all the studied magnetoelectrics except for 
FeTe2C>5Br we introduce praphase structures, which are 
more symmetric crystal structures that can be obtained 
by small atomic displacements. In the cases of MJ1WO4 
and CuO we show that their magnetically ordered phases 
are described by a single IR of the space group of the 
praphase. Therefore, the complex phase diagrams of the 
studied magnetoelectrics can be interpreted as induced 
by a single IR either of the space group of the praphase 
or of the symmetry group of the exchange Hamiltonian, 
which explains the proximity of various magnetic insta- 
bilities in magnetoelectrics. Furthermore, even when the 
introduction of the praphase does not lead to the reduc- 
tion of the number of order parameters, taking into ac- 
count latent pseudosymmetries allows determining the 
influence of external magnetic or electric fields or elas- 
tic stresses of certain symmetry on the phase transitions. 
For MnWC>4 and CuO we obtain temperature-magnetic 
field phase diagrams for certain field directions. The 
phase diagram for MnW04 is in excellent qualitative cor- 
respondence with the experimental one, whereas that for 
CuO predicts new magnetic field-induced phases. The 
magnetic structures of the new phases in CuO as well 
as those of the phases HF, IV, and V of MnW04 are 
discussed. 

The paper is organized as follows. In Sec. |TT] we recall 
the basics of the symmetry of the exchange Hamiltonian 
and its relation to the crystal symmetry, in Sec. IIIII we 
briefly discuss the application of the praphase concept 
to the theory of phase transitions, in Sec. IIVI we per- 
form the analysis of magnetic phase transitions in several 
magnetoelectrics using the exchange symmetry and the 
praphase concept, and in Sec. [V] we discuss the obtained 
results. 



II. EXCHANGE SYMMETRY 

Exchange interaction prevails in the magnetic struc- 
ture formation of most of the magnetically ordered 
crystalsi 30 i 31 In the simplest case the exchange Hamil- 
tonian can be written as 



H 



/ Jnimj 



\&ni ' &mj ) j 



(i) 



where S„i is the ith atom spin operator in the nth unit 
cell and J n imj is the exchange integral between the atoms 
ni and mj. In the system of spins exceeding 1/2 the 
Hamiltonian Eq. (fTJ contains terms of higher order with 
respect to (S n i ■ S mj -). Many-spin exchange is also possi- 
ble. 



The exchange Hamiltonian Eq. ([T]) depends only on 
the mutual orientation of spins and is, therefore, invari- 
ant with respect to simultaneous arbitrary angle rotation 
of all spins around any axis. The orientation of the mag- 
netic structure relative to the crystal axes is determined 
by the relativistic interactions (e. g. spin-orbit or dipolc- 
dipole interaction), which in many cases are weaker than 
the exchange interaction; 30 ! 31 

When building a phenomenological theory of magnetic 
phase transitions in crystals one usually proceeds from 
the symmetry of the disordered paramagnetic phase. The 
experimentally observed magnetically ordered states are 
then classified according to IR's of space group of the 
paramagnetic phase. Thus, the exact symmetry of the 
exchange as well as relativistic interactions is taken into 
account. However, the symmetry of exchange interac- 
tions is higher, than that of the relativistic ones. Con- 
sequently, in cases when the exchange interaction pre- 
vails in determining the magnetic structure, the infor- 
mation on the exchange interaction symmetry is lost. 
Additional symmetry of the exchange interactions leads 
to additional exchange energy degeneracy of magnetic 
states, compared to that corresponding to the space 
group symmetry^ Therefore, one first has to consider 
the symmetry of the exchange Hamiltonian and then take 
into account the exact space symmetry of the crystal. 

The symmetry group of the exchange Hamiltonian 
(the exchange group) is given by the direct product 
G a <8> O s <E> Is, where G a is the space group acting on 
the atom coordinates, O s is the rotational group in the 
spin space, and I s is the inversion group containing the 
unit element and spin inversion.— It can be shown that 
every IR d^ v of the crystal space group, which is char- 
acterized by the wave vector star {k} and which enters 
into the permutational representation of the atoms form- 
ing the magnetic structure, induces IR d^ v ® V of the 
exchange group . 31 ' 32 Here V is the representation giving 
the transformational properties of a pseudovector. To 
determine the connection between the irreducible repre- 
sentation of the exchange group with those of the space 
group one has to expand the limitation of dS^ v <S> V on 
the space group with respect to its IR's^ 2 - 



rfW^l/'.^r^. 



(2) 



Here r" are the expansion coefficients. The irreducible 
representations entering into expansion Eq. ^ are de- 
generate with respect to exchange energy and form the 
exchange multiplet. 32 The splitting of the energy of states 
corresponding to different IR's in the exchange multiplet 
is determined by anisotropic relativistic interactions and 
in many cases is small due to the smallness of these inter- 
actions compared to the exchange energy. Therefore, the 
magnetic structure may be determined by a set of IR's of 
the space group, which form the exchange multiplet /2£r— 



III. THE PRAPHASE CONCEPT 

The praphase concept is widely used in many fields 
of science from elementary particle physics to the solid 
state. In the phase transition theory in crystals the 
praphase concept is widely used, for example, in descrip- 
tion of reconstructive phase transitions ; 33 Adapting it to 
the needs of the present work the concept can be briefly 
described as follows. When building a phenomenological 
theory of phase transitions one starts with a symmetry 
G of the parent phase. Sometimes the crystal structure 
of the parent phase can be transformed by small relative 
displacements of the constituting atoms to a structure 
of higher symmetry G p . One may expect that this high 
symmetry structure, referred to as praphase, is realized 
upon temperature increase, however the decomposition 
or melting of the sample may occur earlier. The phase 
transition G p — > G is described by a generally multicom- 
ponent order parameter {77^}. Thus, one can base the 
phase transition model on the higher symmetry praphase 
structure taking into account that r\i 7^ 0. 

The introduction of praphase improves the phase tran- 
sition model by taking into account latent pseudosymme- 
tries in the object under study. In the following we ana- 
lyze a set of magnetoelectrics using the praphase concept 
and exchange symmetry. The praphase structures of the 
pyroxenes and Cu3Nb2 0s were found using the program 
PSEUDO located at the Bilbao Crystallographic Server 
(http://www.cryst.ehu.es/])> 



34 



IV. PHASE TRANSITIONS IN 
MAGNETOELECTRICS 

A. MnW0 4 

Wolframite M11WO4 has become one of the prominent 
examples of magnetoelectrics since it shows incommen- 
surate paraelectric and ferroelectric magnetically ordered 
phases, low-temperature commensurate magnetic phase, 
and complex temperature - magnetic field phase dia- 
grams. At room temperature it possesses a monoclinic 
structure described by the space group P2/c (C^)* 2 ^ On 
lowering the temperature MnW04 undergoes a sequence 
of magnetic phase transitions at Tn—13.5 K, T2=12.7 K, 
and Ti=7.6 K, which lead to the appearance of magnet- 
ically ordered states AF3, AF2, and AF1, respectively^ 
The structure of the low-temperature commensurate 
magnetic phase AF1 is described by the wave vector 
k c = (1/4; 1/2; 1/2), whereas the incommensurate phases 
AF2 and AF3 are characterized by the wave vector 
k inc = (-0.214; 1/2; 0.457)^In the phases AF1 and AF3 
the collinearly aligned magnetic moments are confined to 
the ac plane forming an angle of about 35° with the a axis 
(this direction is hereafter referred to as the easy axis), 
whereas in the AF2 phase there appears an additional 
component along the b axis. 



Electric polarization along the 6 axis appears at the 
second phase transition at T 2 from the paraelectric AF3 
phase to the ferroelectric AF2 one . 35 ' 36 The polarization 
Pb continuously changes through T2 and drops abruptly 
to zero at T\. The dielectric constant ej, shows a sharp 
peak at T2 and a steplike change at T\ when the ferro- 
electric order disappears ! 35 ' 36 

The magnetic phase transitions in wolframite were 
studied theoretically in a number of works . 16 ' 18 ' 19 ' 37 In 
Refs. [16lfl8l] the authors developed phenomenological 
models of phase transitions in MnWCU starting with 
the order parameters belonging to the incommensurate 
wave vector ki nc . In Ref. [37| only the spin components 
along the easy and b axes are considered, which pre- 
cludes the description of the magnetic field-induced HF 
phase as discussed below i^ 8 . In contrast, in our previous 
work 19 we developed a phenomenological model of phase 
transitions using magnetic order parameters belonging 
to the k c point of the Brillouin zone and accounted for 
the incommensurate phases by means of Lifshitz invari- 
ants. We also suggested that the phase transitions in 
wolframite be described starting from the orthorhombic 
praphasei 19 ' 20 In this section using the praphase concept 
and exchange symmetry we discuss the magnetic struc- 
tures of the magnetic field-induced phases HF, IV, and 
V and the magnetic field-induced flop of electric polar- 
ization (with magnetic field applied along the easy axis). 

In Ref. [19j we developed a phenomenological model 
of magnetic phase transitions in wolframite starting from 
the monoclinic structure P2/c. The magnetically ordered 
phases are described by the order parameters belonging 
to the star of wave vector k c . In this point of the Brillouin 
zone the space group possesses two two-dimensional IR's 
G\ and G^. The phases AF3 and AF1 are described by 
Gi only, whereas G\ additionally condenses in the phase 
AF2.2 5 - (It has to be noted that in Refs. [l^-[2l[ we have 
inadvertently assumed that the phases AF3 and AF1 are 
described by G\ only, whereas Gi additionally condenses 
in the phase AF2. However, this does not change any of 
the obtained results.) 

Therefore, upon lowering the temperature MnW04 ex- 
periences two close magnetic instabilities with respect 
to different IR's G\ and G2. It was argued that the 
phase transitions in wolframite can be described us- 
ing orthorhombic praphasei 19 ' 20 In the monoclinic struc- 
ture with the monoclinic angle f3 s=s 91° the atoms 
are located at positions Mn - (1/2,0.6853,1/4), W - 
(1/2,0.3147,3/4), Oi - (0.2108,0.1024,0.9419), and 2 
- (0.2516, 0.3752, 0.3931). 25 The displacement of oxygens 
towards positions of higher symmetry Oi - (0.2108,0,0) 
and 2 - (0.2516, 1/2, 1/2) and setting (3 = 90° results in 
the orthorhombic structure described by the space group 
Pmcm (-D 5 ^). I n t ne following we define the orthogonal 
coordinate axes x, y, and z parallel to the crystal axes 
a, b, and c of the orthorhombic praphase, respectively. 
The phase transition Pmcm-P2/c is described by the 
homogeneous deformation tensor component U xz , which 
should be taken nonzero when building a model of phase 



transitions in MnWC>4. 

The modulation wave vector k c retains its position in 
the orthorhombic structure. In this point of the Bril- 
louin zone the space group Pmcm possesses one four- 
dimensional IR Pi. Thus, the IR's G\ and Gi stem 
from Pi, which splits under the action of U xz . Therefore, 
the magnetically ordered states of MnW04 can be inter- 
preted as induced by a single IR P\ of the orthorhombic 
praphasei 19 ' 20 

The magnetic representation analysis starting from the 
praphase can be performed as follows. The magnetic mo- 
ments of the two Mn 2+ ions in MnW04 can be expressed 
in the form 



Mi = 



M, 









e lklt + 




Ak 9 -.L 
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~ikit 



where t is the lattice vector, k\ — k c , and &2 = —k\. The 

k ■ 

four quantities M iQ | for every direction a transform ac- 
cording to the IR Pi, i.e., Pi enters three times into the 
magnetic representation of Mn 2+ ions. Thus, three order 
parameters (510,32a, 53a, ff4a) (a = X,y,z), which give 
the magnetic moment components along a and trans- 
form according to Pi, determine the magnetic structure 
in wolframite. 

At low temperatures in the phase AF1 (and also in the 
phase HF as suggested in Ref. [38]) the magnetic mo- 
ments are confined to the ac plane. Therefore, in order 
to determine the magnetic structure of the HF phase we 
consider the order parameters with a = x and z only. 
The quadratic part of the thermodynamic potential ex- 
pansion with respect to these order parameters can be 
written in the form 



F 



-{Ix+Iz)+fIxz + Kl L 



(3) 



where a, /, and k are phenomenological coefficients, I x = 

9ix + six + si* + sL> h = g\ z + sL + si + g\ z , hz = 
gixgiz + g2 X g2z + g3xg3z + g4xg4z, and l v = U xz (gi x gi z + 

gixg2z - gzxgsz - g^xgAz)- The term proportional to n, 
reflects the splitting of the order parameters due to the 
monoclinic distortion (U xz ) and external magnetic field 
H applied in the ac plane. Since H X H Z transforms as U xz 
under the symmetry elements of the space group and, 
therefore, directly influences the splitting of the order 
parameters, we assume U xz — U xz — wH x H Zl where w 
is a phenomenological coefficient. In Eq. ^ we use the 
same coefficient at I x and I Z1 since the order parameters 
gi x and gi z belong to the same exchange multiplet^L The 
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FIG. 1: (a) and (d) give two basic ordering patterns of Mn 2+ 
spins corresponding to the exchange multiplets G\ ® V and 
G2 <S> V, respectively. The plus and minus signs give relative 
spin directions, (b) and (c) show magnetic structures in the ac 
plane with different directions of the easy axis corresponding 
to the pattern shown in (a), whereas (e) and (f) correspond to 
(d). The magnetic structures (b), (c), (e), and (f) ofMnWOi 
correspond to qi 7^ 0, 34 7^ 0, qa 7^ 0, and q$ 7^ 0, respectively. 



diagonalization of the quadratic form Eq. ([3]) results in 

F = (a + f-KU' xz )(ql + q 2 2 ) 
+(a-f + KU' xz )(q 2 3 +ql) 
+(a-f-KU' xz )(ql+ql) 

+ (a + f + K U xz )(q 2 7 + q 2 a ) 1 (4) 

where qi = (g ix + gi z )/2, 92 = (33a; + #3z)/2, qz = 
(g4z-g4x)/2, qi = (g 3 z-g3x)/2, q 5 = {g2 Z -g2x)/2, q 6 = 
(giz-9ix)/2, 17 = (92x + g2z)/2, and q 8 = {gi x +gi z )/2. 
The experimentally observed magnetic structure of the 
phase AF1 is shown in Fig. QJb).— The easy axis, along 
which the magnetic moments are aligned in the AF1 
phase, makes an angle of 35° with the a axis. It has 
to be noted that in Eq. ([3]) we for simplicity did not 



include the invariants U' xz (g\ x + g\ x — g 2 x — g 2 x ) and 
U'xzidlz + dlz - glz - ff£J> which determine the angle 
between the easy axis and the a axis. The magnetic 
structure of the AF1 phase (i.e., at zero and low fields) is 
described by the solution q 2 7^ 0. Therefore, from Eq. Q 
it can be argued that / < and kU xz > 0, which makes 
this state preferable at zero applied field and T < T\. 
(Here we assume that n > and U xz > without loss of 
generality.) 

The experimental data reveal that below T\ the ap- 
plication of the magnetic field along the easy axis re- 
sults in the field-induced AF1-AF2, AF2-HF, HF-IV, 
IV-V, and V-PM sequence of phase transitions^Sr— In- 
formation on the magnetic structure of the HF phase 



is scarce. However, neutron diffraction studies suggest 
that the field-induced HF phase is commensurate with 
the same modulation vector as in the phase AF L 38 ' ° 
It is argued that HF is a spin-flop phase, in which the 
magnetic moments are switched perpendicular to the ap- 
plied field within the ac plane with the same relative spin 
arrangement as in AFl^ A simple uniform rotation of 
the spins would result in magnetic structure shown in 
Fig. [TJc), which is given by (74 7^ 0. In our model we 
assume that the effective monoclinic splitting U xz is di- 
rectly influenced by the magnetic field applied in the ac 
plane such that H X H Z 7^ 0. Therefore, it can be argued 
that sufficiently high magnetic field parallel to the easy 
axis (i.e., H X H Z > 0) changes the sign of the splitting 
U xz (thus, w > 0). Consequently, assuming that loH x H z 
is the main contribution to the thermodynamic poten- 
tial, the solution q$ 7^ will become stable instead of 
(72 7^ as follows from Eq. (U|). The corresponding mag- 
netic structure is shown in Fig. [TJf ) and is different from 
the one shown in Fig. [IJc). 

Using the monoclinic lattice as the reference point the 
magnetic structures of wolframite can be described by 
two exchange multiplets generated by the IR's G\ and 
G2 2^ The respective spin ordering patterns are shown in 
Figs -da) andQJd) an d can not be transformed into each 
other by the monoclinic symmetry group operations. The 
magnetic structure of the phase AF1 corresponds to the 
ordering pattern of Fig. [TJa). This pattern is given by 
the exchange multiplet Gi<E>V, whose limitation on the 
space group splits into G\ © 2G2 , which reflects the fact 
that both x and z components of magnetic moments are 
described by IR G2, whereas the y component by G±. 

The magnetic structure of the phase HF [Fig. [Tic)], 
which is suggested in Ref. [38(], corresponds to the same 
exchange multiplet as that of AF1. However, using the 
orthorhombic praphase as reference we find that the two 
multiplets G\ ® V and G2 <8> V stem from a single mul- 
tiplet Pi <g> V generated by the IR Pi of the praphase. 21 
The splitting of the exchange energies of the magnetically 
ordered states corresponding to multiplets G\ <g> V and 
G2 <£) V is determined by effective monoclinic distortion 
U xz , which is directly influenced by the magnetic field 
with H X H Z 7^ 0. Therefore, sufficiently high magnetic 
field along the easy axis changes the sign of U xz , which 
results in lower exchange energies of magnetically ordered 
states corresponding to the exchange multiplet G2 <8> V . 
Therefore, the magnetic structures of the phases AF1 
and HF are described by different exchange multiplets 
{G\ <g> V and G2 ® V, respectively). Thus, we argue 
that the phase HF possesses magnetic structure shown 
in Fig.QJf), which is described by IR G\. 

The fact that U' xz changes sign at magnetic field H xz sa 
8 — 10 T along the easy axis is supported by shrinkage of 
the temperature range of stability of the phase AF3 at 
H xz , which results in direct transition from the param- 
agnetic phase to the phase AF2 i 40 ' 41 Therefore, U xz = 
at H xz , which implies the absence of splitting of IR Pi 
allowing for simultaneous condensation of IR's G± and 



G 2 directly from the paramagnetic phase to the phase 
AF2. At magnetic fields H xz > H xz the splitting U xz 
has sign different from that at H xz < H xz , which should 
favor the condensation of the IR Gi first upon lower- 
ing the temperature from the paramagnetic phase. This 
implies that the magnetic order in the experimentally 
observe d 40 ! 41 phase V is described by IR Gi, in contrast 
to G 2 for the phase AF3. Further temperature lowering 
at H xz > H xz should result in additional condensation of 
Gi and a phase transition to either AF2 or IV depending 
on H xz . 

According to neutron diffraction experiments, in the 
phase AF3 sinusoidally modulated magnetic moments lie 
in the ac plane forming an angle with the a axis similar 
to that in the commensurate AFl phase. 25 Similarly, the 
preceding analysis argues that the phase V is incommen- 
surately modulated with magnetic moments lying in the 
ac plane perpendicular to their direction in the phase 
AF3 and ordered according to the exchange multiplct 
Gi®V [Fig. UK)]. 

To support this interpretation we perform numeri- 
cal minimization of the following thermodynamic po- 
tential expansion with respect to the order parameter 
(ci, ci, c 3 , C4) transforming according to IR Pi 
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where the invariants are I\ = c\ -f- cf + c; 



h = 



+ ct + c 



-1 T L 2 t u 3 

ifi = C4C3 + C2C4 
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C1C2C3C4, 



h 



3 1" ^4! 

Ci Cn "T" CqC^ 



1^2 t K3°4' 
5 = CJCg I" C2CJ, /„ = U' xz {c 1 C 3 + C 2 C 4 ), I p = P y (ciC 2 - 

C3C4), Ip = H x H z U xz P v (cic 2 - C3C4), I L = cidci/dx - 
Cidci/dx + c-idcijdx — cidc^/dx, and Is — (dci/dx) 2 + 
(dc 2 /dx) 2 + (dc 3 /dx) 2 + (dc 4 /dx) 2 , and V is the volume 
of the sample. In the thermodynamic potential Eq. (|5|) , 
which is invariant with respect to the space group ele- 
ments of the praphase, we for simplicity consider the de- 
pendence of the order parameter on x only. The following 
values of phenomenological coefficients v = 15, B\ = 20, 
B 2 = 4, B 3 = -1, B 4 = -21, B 5 = -1, w = 1, a = 1, 
(5 = 1, U xz = 0.09, k = —40 and a=500 give the phase 
diagram shown in Fig. Ufa). In the Landau theory of 
phase transitions it is usually assumed that the coefficient 
at Ii possesses the strongest dependence on external pa- 
rameters such as temperature or external magnetic field. 
Therefore, the A axis in the phase diagram of Fig. HJa) 
can be associated with temperature. We also included 
the term proportional to v, which describes the simplest 
magnetic field influence and results in the decrease of 
temperatures of phase transitions associated with IR P\ . 
It has to be noted that the phenomenological constant a 
is related to the dielectric susceptibility % by a = x~ l ■ 

Performing the substitution of variables according to 
ci = -gi + 92 + g 3 + .94, c 2 = .91 + g 2 - g 3 + g±, c 3 = 
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FIG. 2: (color online) (a) Calculated phase diagram of 
MnWC>4 for magnetic field applied along the easy axis. 
Dashed and solid lines represent phase transitions of the sec- 
ond and first kind, respectively. Dotted line represents the line 
at which the electric polarization changes sign, (b) Electric 
polarization dependence on the magnetic field H at constant 
A. The red, green, and blue lines correspond to the respective 
paths shown in (a) by colored axes. 
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2 — g\), which reflects the splitting of the 



-54, and C4 

2C/L(5I+52 2 -53 2 
order parameter (51,52,53,54) into two different order 
parameters (51,52) and (53,54) under the influence of 
U' xz . It can be shown that (gi,g 2 ) and (53,54) transform 
according to IR's G\ and G 2 of the space group P2/c, 
respectively. The magnetoelectric interaction takes the 
form I p = 4P y (gig 3 + 3254), which reflects the necessity 
of condensation of both IR's G\ and G 2 in order for P y 
to arise. 

The phase diagram Fig. |^a) obtained by minimizing 
the functional Eq. (O is in excellent qualitative agree- 
ment with the experimentally observed diagram for mag- 
netic field applied along the easy axis42r— At zero mag- 
netic field M11WO4 shows the sequence of phase transi- 



tions AF3-AF2-AF1 at decreasing temperature. Appli- 
cation of magnetic field along the easy axis results in 
shrinkage of the temperature range of stability of the 
phase AF3. At U' xz = 0, which occurs at H c xz = 0.3 
for the taken values of phenomenological coefficients, the 
phase AF3 disappears corresponding to the absence of ef- 
fective monoclinic splitting. This value of magnetic field 
corresponds to the widest temperature range of stability 
of the phase AF2. Further magnetic field increase results 
in the growth of U xz in absolute value but with differ- 
ent sign, which leads to the appearance of the phases V 
and HF, the structures of which were discussed above. 
The magnetic phases in the phase diagram Fig. HJa) 
are given by the following values of the order param- 
eter (51,52,53,54): AF3 - (0,0,g 3 (x),gi(x)), AF2 and 
IV - (gi(x),g 2 (x),g 3 (x),g4(x)), AF1 - (0, 0,53,0), V - 
(51 (x), 52(z),0,0), and HF- (51, 0,0,0). 

The experiment reveals the field-induced sign change 
of electric polarization Pb for magnetic fields along the 
easy axis42 which has led to the introduction of the new 
phase IV.— ~— From our point of view such behavior can 
be explained by the sign change of the coefficient at the 
magnetoelectric invariant. In Eq. ([5]) we take this into 
account by the invariant /' , which together with I p deter- 
mines the magnetoelectric response. The resulting mag- 
netic field dependencies of Pb at different values of A are 
given in Fig. HJb) and are in good qualitative correspon- 
dence with the experimental ones42 

The experimentally observed wave vector of the incom- 
mensurate phase IV (—0.215, 0.503, 0.460) is close to that 
of AF24S The appearing modulation along the b axis 
can be accounted for by the invariant H x H z {c\dc±/ dy — 
Cidci/dy + C2dc 3 /dy — c^dc2/dy), which is proportional 
to the magnetic field with H X H Z 7^ 0. Therefore, the 
magnetic field with H X H Z ^ also induces spatial mod- 
ulation of the magnetic order along the y axis due to the 
Lifshitz invariant proportional to H X H Z . 

The experiment reveals the memory effect in MnWC^ 
for the magnetic field-induced phase transitions sequence 
AF2-HF-IV, which consists in the fact that irrespective 
of the poling electric field direction the directions of the 
electric polarization Pb in the ferroelectric phases AF2 
and IV are always opposite to each other42 This ef- 
fect persists also when the thermodynamic path crosses 
the non- ferroelectric phase HF. According to our model 
the phase HF is paraelectric and the information on the 
previous polarization direction should be completely lost 
upon the phase transitions AF2-HF or IV-HF. Therefore, 
considering the field-increasing run, the subsequent po- 
larization direction in the phase IV when reached from 
the phase HF should be determined by the applied elec- 
tric field, which is not the case in the experiment. How- 
ever, from our point of view this memory effect is dy- 
namic. Indeed, it is found that the characteristic time of 
the response of the magnetic order parameter in MnWCU 
to the electric field pulses is of the order of 10 ms4£ 
This time span is of the order (or even longer) then that 
of the magnetic field pulses used to study the memory 



effect^ Therefore, the interpretation, which takes into 
account possible preservation of embryos of the ferroelec- 
tric phased together with the sign change of the magne- 
toelectric coefficient suggested in our model may explain 
this memory effect. Experiments with longer magnetic 
field pulses are required in order to confirm or refute the 
dynamic nature of the memory effect. 

Thus, we have built a model describing the behavior 
of MnWC>4 in magnetic field applied along the easy axis. 
Using the praphase concept we identified the influence 
of the magnetic field applied parallel to the easy axis on 
the thermodynamic potential expansion, which consists 
not in simple quantitative field dependence of the phe- 
nomenological coefficients, but in new terms allowed by 
the symmetry. This approach also allowed us to suggest 
possible magnetic structures for the phases HF and V. 

It has to be noted, that the phase diagram of Fig. HJa) 
is topologically similar to that published earlier^ In the 
present work, however, we slightly modified our model 
and adjusted the phenomenological constants in order to 
account for new details such as the suppression of the 
antifcrromagnetic order by magnetic field, and magnet- 
ically induced flop of electric polarization and the AF1- 
AF2-HF-IV-V phase transition sequence. 

In order to study the problem qualitatively we for sim- 
plicity carried out the expansion of the thermodynamic 
potential Eq. ([5]) only up to the fourth order with re- 
spect to the magnetic order parameter and included only 
the required magnetic field dependence. The experimen- 
tal phase diagrams of wolframite are characterized by 
the values (Tjy — T)/Tn close to 1. Therefore, better 
quantitative comparison with the experiment requires ex- 
pansion of the thermodynamic potential to higher orders 
and taking better account of the magnetic field influence, 
which is beyond the scope of this paper. 



B. CuO 

Cupric oxide CuO is a type-II magnetoelectric with 
one of the highest Curie temperatures— At normal con- 
ditions it possesses a monoclinic structure described by 
the space group C2/c (C% h )£^ Upon lowering the tem- 
perature CuO undergoes three magnetic phase transi- 
tions at Tjv 3 =230 K, T N2 =229.3 K, and T N1 =213 K, 
which lead to the appearance of magnetically ordered 
phases AF3, AF2, and AF1, respectively^ The phase 
AF2 is characterized by long-wavelength modulation 
with the wave vector (0.506, 0, —0.483), whereas the com- 
mensurate phase AF1 is described by the wave vector 
(1/2,0,-1/2)4241 The existence of the phase AF3 was 
confirmed only recently by ultrasonic measurements and 
no experimental data on its magnetic structure exist to 
date4& Electric polarization in CuO along the b axis ap- 
pears in the incommensurate phase AF2.— The newly 
discovered phase AF3 is argued to be paraelectric since 
no anomaly of the dielectric constant is observed at Tjv3 
in contrast to T/V2- 46 
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Magnetic phase transitions and magnetoelectricity in 
CuO were studied theoretically in a number of works. 
Monte-Carlo studies based on first-principles calculations 
failed to reproduce the existence of the phase AF3 , 48 ' 49 A 
phenomenological theory of phase transitions in CuO was 
suggested in Ref. [26| , in which the phase transition PM- 
AF2 was described using the triggering mechanism. This 
implies the first order character of the PM-AF2 transi- 
tion. In contrast, earlier based on a phenomenological 
theory we suggested that similar to M11WO4 there should 
be an intermediate phase between the paramagnetic and 
ferroelectric phases^ We also suggested an orthorhom- 
bic praphase for the description of CuO. Here we further 
develop this approach to CuO suggesting possible mag- 
netic structure of the phase AF3 and building the phase 
diagram for magnetic helds applied in the ac plane at an 
angle to both crystal axes, which shows two new field- 
induced phases. 

In the monoclinic structure (/3 « 99°) the copper 
ions occupy position (i, j, 0) and oxygen ions - (0, y, i) 
(y = 0.416) 4£ It can be found that y = ± and /3 = 90° re- 
sult in the orthorhombic structure described by the space 
group Cccm (Df/l)- 20 Similar to wolframite the phase 
transition Cccm-C2/c is described by the deformation 
tensor component U xz , which transforms according to 
IR GM 4+ of the symmetry group of the praphase and 
has to be assigned nonzero value in order to describe the 
monoclinic phase. 

The low-temperature magnetic structure AF1 is de- 
scribed by the wave vector k c = (4,0,-4), whereas 
the incommensurate phase AF2 by (0.506,0,-0.483), 
which is close to the commensurate valuei 44 ' 47 There- 
fore, one can describe the phase transitions in CuO using 
the commensurate wave vector k c and account for long- 
wavelength modulation by Lifshitz invariants allowed by 
the symmetry. 20 Using the monoclinic structure as ref- 
erence it was shown that the magnetic structure of the 
phase AF2 is described by two IR's Bi and B2, whereas 
the phase AF1 by B 2 m±2£. 

The magnetic representation analysis using the or- 
thorhombic praphase can be performed as follows. In 
the following we define the orthogonal coordinate axes x, 
y, and z parallel to the crystal axes a, b, and c of the 
orthorhombic structure, respectively. The wave vector 
k c maintains its position in the orthorhombic structure. 
Thus, the magnetic moments M\ and M 2 of the two cop- 
per atoms Cui and Cu2 present in the primitive cell and 
located at positions (|,-|,0) and (|, §, 3), respectively, 
can be expressed as 



M„ 



M kl 

nx 

M kl 

ny 

M kl 
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M k2 
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M^ 2 

±V - L ny 
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where n = 1 or 2, t is the lattice vector, k\ — k c , and 

7* 7* k ■ 

ki = —k\. The four quantities M n 3 a for every direction a 
transform according to IR A± of the space group Cccm, 



i.e., A\ enters three times into the magnetic represen- 
tation of Cu 2+ ions. This implies that the two two- 
dimensional IR's Bi and B 2 of the monoclinic space group 
merge into a single four-dimensional IR A\ of the space 
group of the praphase. It explains the closeness of insta- 
bilities with respect to Bi and B2 on the thermodynamic 
path and allows interpreting the magnetically ordered 
states of CuO as induced by a single IR Aii 20 ' 21 Thus, 
the magnetic representation analysis for CuO closely re- 
sembles that for M11WO4 preformed in Sec. IIV Al In fact 
the analogy can be drawn much further. 

Let us denote by (31,32,33734) the magnetic order pa- 
rameter transforming according to A\. It can be shown 
that (171,32) and (33,174) transform according to Bi and 
B2, respectively, when only the elements of the mono- 
clinic space group are considered. Similar to MnW04 
this is reflected in the existence of the invariant 

KM + al- oi-al), 

which splits the order parameter gi^ Here we again use 
U xz = U xz — wH x H z since H X H Z transforms as U xz un- 
der the symmetry elements of the praphase. The mag- 
netoelectric interaction responsible for the appearance of 
electric polarization P y in the phase AF2 is given by 



I p = P y U x , z (g 1 g3 + g 2 g4), 



(0) 



where U xz = U xz — w p H x H z and w p is a phenomenologi- 
cal constant. Therefore, condensation of both B\ and B^ 
is necessarily for its appearance. As follows from Eq. (|6]) 
the electric polarization P y oc {U xz — w p H x H z ) is directly 
influenced by the magnetic field with H X H Z 7^ 0. Thus, 
depending on the sign of H x H Zl which can be changed by 
choosing suitable field direction, one may expect different 
magnetic field dependence of P y . For example, if P y > 
then one may either expect its further increase with the 
magnetic field or, similar to MnW04, its decrease with 
subsequent sign change (i.e., a polarization flop). 

The thermodynamic potential expansion with respect 
to the order parameter (ci, €2,03, C4), which transforms 
according to A\ and is related to gt (i=l,2,3,4) by c\ = 
-.91+32+33+34, c 2 = .91 +.92 -.93 +34, c 3 = -31+.92-33- 
34, and C4 = #1 + 32 + .93 — 34, can be written in the form 
of Eq. ([5]) with the difference that I p is given by Eq. (J5]) 
and I' = 0. The following values of phenomenological co- 
efficients v = 4, Bx = 20, B 2 = 4, B 3 = -1, B 4 = -21, 
B 5 = -1, w = 1, w p = 0, a = 1, 6 = 1, U xz = 0.14, 
and a=500 give the phase diagrams shown in Fig. [3] 
The order parameter takes the following values in the 
magnetically ordered phases: AF3 - (0, 0,(73(2), 34(2;)), 
AF2 - ( gi (x),g 2 (x),g 3 (x),g 4 (x)), AF1 - (0,0,0, g), HF - 
(0,3,0,0), and AF3' - (gi(x), g 2 (x), 0,0). We find that 
the phases AF3, AF2, and AF3' are incommensurate 
and are characterized by the same modulation vector, 
whereas the phase AF2 is also ferroelectric with P y 7^ 0. 
Similar to wolframite the phases AF3 and AF1 are de- 
scribed by IR B2, whereas HF and AF3', which corre- 
sponds to the phase V of MnW04, by B\. The phase 
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FIG. 3: (a) and (b) give calculated magnetic phase diagrams 
of CuO for magnetic fields parallel to [101] and [101], i.e., for 
H X H Z > and H X H Z < 0, respectively. Dashed and solid 
lines represent phase transitions of the second and first kind, 
respectively. 



diagram of Fig. (3Ja) resembles that of M11WO4 shown in 
Fig. [2ja) due to the similarity in the thermodynamic po- 
tential expansion. The external magnetic field along [101] 
changes the splitting of the order parameter gi , which for 
sufficiently high fields results in appearance of the phases 
AF3' and HF, which are described by IR B\. In contrast, 
the magnetic field parallel to [101] strengthens the split- 
ting resulting in shrinkage of the temperature range of 
stability of the phase AF2, which is described by both 
B\ and B2, and widening of the stability ranges of the 
phases AFl and AF3, which are described by B^. The 
phase diagrams of Fig. [3] are obtained under the assump- 
tion w > 0. If w < 0, then Figs.[3ja) andEJb) correspond 
to the magnetic fields applied along [101] and [101], re- 
spectively. 

The magnetic representation analysis similar to that 
performed for MnWC>4 in Sec. IIV Al allows determining 
the magnetic structures appearing in the phase diagrams 
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FIG. 4: (color online) (a), (b), (c), and (d) give the calcu- 
lated magnetic structures AFl, AF3, HF, and AF3' of CuO, 
respectively. The magnetic moments are directed along the b 
axis. 



of Fig. [3] Figure 2] presents the magnetic structures AFl, 
AF3, HF, and AF3'. In the phase AF3 the Cu 2+ mag- 
netic moments are sinusoidally modulated and directed 
along the b axis. It has to be noted, that according to 
our model all of the magnetic moments order in the phase 
AF3 as follows from Fig.QJb) in contrast to the result ob- 
tained in Ref . [46| . From our point of view the authors of 
Ref. [46| find that only half of the moments order in the 
phase AF3 due to the fact that in their model the mag- 
netic structures of CuO are described by complex pseu- 
dovector S. This implies that for every direction a = x, 
y, and z the magnetic structure is described by two or- 
der parameters ReS^ and Im5 Q , whereas as we showed 
above the magnetic representation for Cu 2+ ions for ev- 
ery direction a is four-dimensional. 

Thus, our model of phase transitions in CuO predicts 
the magnetic field-induced phase transition AF1-HF for 
#|| [101] (or for #|| [101] if w < 0). In MnW0 4 in a sim- 
ilar phase transition (although through the intermediate 
phase AF2) the easy axis changes from being parallel 
to the applied field to perpendicular direction. There- 
fore, one is tempted to interpret this phase transition as 
a simple spin-flop phase transition, at which the spins 
prefer to align along the direction perpendicular to the 
applied magnetic field. However, as we showed above in 
Sec. IIV Al the mechanism of this phase transition can be 
more complex and deserves additional study of the micro- 
scopic causes of the differences in the exchange energies 
corresponding to two exchange multiplets. In contrast, 
in CuO as shown in Figs.HJa) and0Jc) this phase tran- 
sition is between two phases, for both of which the easy 
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axis is the b axis, i.e., perpendicular to the applied mag- 
netic field. Therefore, experimental confirmation of the 
suggested magnetic field-induced phase transition will 
strongly support the suggested models of the phase tran- 
sitions in both MnWC>4 and CuO. 



C. Pyroxenes NaFeSi 2 O e and NaFeGe 2 O e 

Three members of the pyroxene family (with the 
general formula AMSi 2 6 ) NaFeSi 2 6 , LiFeSi 2 6 and 
LiCrSi 2 06 were recently shown to display ferroelectric- 
ity induced by magnetic order. 50 In this work we fo- 
cus on NaFeSi 2 6 (NFSO) and NaFeGe 2 6 (NFGO). In 
the paramagnetic phase both compounds possess mon- 
oclinic symmetry described by the space group C2/c 
(C|fr)] 51 i 52 Upon cooling the compounds show two mag- 
netic phase transitions at Tat = 8 K and T c = 6 K 
for NFSO 50 and T N = 13 K and T c = 11.5 K for 
NFGO 52 leading to the magnetically ordered phases AF2 
and AF1, respectively. Below T/v the magnetic structure 
is found to be incommensurate with k = (0, 0.77, 0)^ and 
k = (0.3357, 0, 0.0814) 52 for NFSO and NFGO, respec- 
tively. NFSO is found to be ferroelectric below T c with 
polarization P||6, 50 whereas in NFGO the polarization 
also appears below T c £l 



1. NaFeSi 2 6 

A phenomcnological model of phase transitions in 
NFSO was suggested earlier using the order parameters 
belonging to the wave vector (0,0.77, O)^ In contrast, 
we now build a model assuming the instability in the 
k = (0,1/4,0), which is close to the equivalent vector 
(0, 0.23, 0), point of the Brillouin zone and account for the 
spatial modulation by considering the Lifshitz invariants 
allowed for the respective order parameters. In this point 
of the Brillouin zone the space group C2/c possesses two 
two-dimensional IR's Ai and A 2 . In the monoclinic C2/c 
structure the iron ions Fei and Fe 2 are located in po- 
sitions (0,0.8991,1/4) and (0,0.1009,3/4), respectively. 
The permutational representation for these ions is given 
by 2Ai. Thus, there are two exchange multiplets given 
by Ai (g> V — Ai © 2A 2 . In the following we define the 
orthogonal x, y and z axes parallel to the a axis, parallel 
to the b axis and perpendicular to both the a and b axes 
of the monoclinic cell, respectively. It can be shown that 
the y components of iron spins transform according to 
Ai , whereas the x and z components according to A 2 . In 
the magnetically ordered phases the spins are confined 
to the ac crystal planei 51 ' 53 Thus, the magnetic order is 
described by IR A 2 . In order to account for the electric 
polarization Pb one should consider two order parame- 
ters (r]x,£x) an d (r] z ,£, z ) transforming according to A 2 
and describing the iron spin components along the x and 
z axes, respectively. The magnetoelectric interaction is 



TABLE I: Atomic positions of NaFeSi2C>6 in the monoclinic 
C2/c and orthorhombic Cmcm structures. In the Cmcm 
structure the monoclinic angle ft — 107.3° should be set to 
,5 = 90°. 
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then given by (r] x £z — £ x r] z )Py, whereas a single order pa- 
rameter (rj x ,£ x ) or (t] z ,^ z ) would require invariants of the 
form Vp£p(VB- 7 '>lj(ip + '?'nptj-£p)Py, where ft = x,z and 
which are of 8th order with respect to the magnetic order 
parameter. The long-wavelength modulation along the y 
axis should be accounted for by the Lifshitz invariants 
rjfsd^/dy-^dTjp/dy. 

However, the magnetic phase transitions in NFSO al- 
low another interpretation using the praphase concept. 
Small atomic displacements towards positions of higher 
symmetry result in the orthorhombic crystal structure 
with Cmcm (D^) symmetry^ The atomic positions in 
the monoclinic 55 C2/c and orthorhombic Cmcm struc- 
tures of NFSO are given in Table |U It has to be noted 
that the monoclinic modification of the pyroxene-type 
MnGe0 3 is isostructural to both NFSO and NFGO 
and possesses significantly smaller monoclinic angle ft ra 
101. 5°£§. The phase transition Cmcm-C2/c is described 
by the component U xz of the homogeneous deformation 
tensor, which transforms according to IR GM 4+ of the 
Cmcm space group. Thus, non-zero U xz should be taken 
into account when describing the phase transitions in 
NFSO starting with the orthorhombic praphase. 

In the k = (0, 1/4, 0) point of the Brillouin zone, which 
retains in the orthorhombic structure, the space group 
possesses four two-dimensional IR's Aj (i — 1,2,3,4). 
The permutational representation of the iron ions is given 
by 2Ai. The limitation of the exchange multiplets, which 
are given by Ai <£> V , on the space group splits into the 
direct sum A30A 2 ©A4, with the x, y, and z spin compo- 
nents transforming according to A 3 , A 2 , and A4, respec- 
tively. Thus, one should consider the order parameters 
{O'x^bx), (a y ,b y ), and (a z ,b z ) transforming according to 
three different IR's A3, A 2 , and A4, respectively. The 
fact that the order parameters (rj x , £ x ) and (rj z , £ z ) in the 
monoclinic structure transform according to a single IR 
is reflected in the existence of the invariant 



lu = ia x b z - b x a z )U xz . 

The magnetoelectric interaction is given by 

I P = (a x a z + b x b z )U xz P y , 

which is, thus, proportional to the deformation U xz . 

The expansion of the thermodynamic potential up to 
the fourth order in powers of the order parameters can 
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be written in the form 

* = V I { T Jl + T /2 + Klu + fl1 ' + h % 



-hJ\ 



1J2 + (T\Il\ + 0"2-Tl2 + 5\Is\ + 82IS2 



-sip + ^Py^<iy 



(7) 



where A 1: A 2 , k, f x , f 2 , f 3 , fi, e\, er 2 , <Si, S 2 , s, and 
A are phenomenological coefficients, I\ = a 2 + b 2 , I 2 — 
a\ + b 2 z , Ji = (a x a z + b x b z ) 2 , J 2 = {a x b z - b x a z ) 2 , 7 L i = 
a x db x ldy-b x da x jdy, I L2 = a z db z /dy~b z da z /dy, hi = 
(da x /dy) 2 + (db x /dy) 2 , and I S2 = (da z /dy) 2 + (db z /dy) 2 . 
In order to minimize the functional Eq. (|7|) we in- 
troduce polar coordinates a x — r x coscj), b x — r x sm(j), 
o-z = Tz cos((fi + A<j>) and b z = r z sin(0 + A0) and assume 
that only the phase <f> is spatially dependent. Minimiza- 
tion of Eq. ([7]) shows that the phase transition from the 
paramagnetic phase to the phase with r x ^= 0, r z =/= 0, 
P y = 0, A<p = 7r/2, and cf> = —qy with 



2Si 



(aid 2 - ct 2 5i)51k 2 U, 



■2 

xz 



2(A 2 S( + ai{ai8 2 - 2a 2 5 1 )) 2 



occurs at 



Ai = 



swu 2 z 



61 A 2 S{ + ai(aiS 2 - 2a 2 Si) 



(8) 



where we used the expansion with respect to kU xz up to 
the second order. We associate this phase with the phase 
AF2. 

In the ferroelectric phase AF1 (r x ^ 0, r z ^ 0, P y ^ 
and Acf> ^ ±tt/2) the polarization is given by P y = 
—sU xz r x r z cos(A(f)/A with 



sin A0 



AkU xz 



r x r z (2Af 3 ~2AU-s 2 U 2 z )- 



(9) 



Thus, the phase transition line to this phase is deter- 
mined by the condition |sinA0| = 1, which is satisfied 
when \r x r z \ becomes sufficiently large. The temperature 
range of stability of the AF2 phase is, therefore, propor- 
tional to the distortion U xz . Similar to the above cases 
of MnW04 and CuO the magnetic field in the ac plane 
with H X H Z =/= will introduce additional contribution to 
the monoclinic distortion U xz , which will affect the in- 
variants Ijj and Ip resulting in the field dependence of 
the transition temperatures T N and T c through Eqs. © 
and ((9|), respectively. Thus, for certain field direction 
and value one may expect zero effective monoclinic dis- 
tortion, which will result in shrinkage of the phase AF2 
and a direct transition to AF1 from the paramagnetic 
phase. The external magnetic field applied in the ac 
plane along (b x c) has indeed strong influence on the 
phase transition temperatures TV and T c and leads to 
the suppression of Pt, and appearance of P c £^ In con- 
trast, the magnetic field along b has little impact on the 
magnetic properties and electric polarization. 



The appearance of P c for -ff||(b x c) can be ex- 
plained by switching of the spin rotation plane from ac 
to be and by the existence of magnetoelectric interactions 
(a y b z — b y a z )P z and (a y b z — b y a z )U xz P x . The coefficient 
at the latter in the thermodynamic potential expansion 
is probably very small since no -P(bxc) is observed. 50 



2. NaFeGe 2 6 

NFGO shows incommensurate magnetic structure at 
2.5 K with the modulation vector (k x , 0, k z ), where k x = 
0.323 and k z — 0.08. 55 It is found that upon the appear- 
ance of the magnetic structure at TV its modulation vec- 
tor is given by k x « 0.295 and k z — 0.065. Therefore, on 
lowering the temperature the k x modulation vector com- 
ponent grows towards the commensurate value l/3i 55 i 57 
The commensurate phase, however, is not realized in 
NFGO. Thus, in order to develop the phenomenological 
theory of phase transitions in it we use the k = (1/3, 0, 0) 
wave vector. In this Brillouin zone point the space group 
C2/c has two two-dimensional IR's B\ and B 2 and the 
permutational representation of the ions Fei and Fe 2 is 
given by B\@B 2 . The magnetic structure in NFGO cor- 
responds to the B 2 <g> V' exchange multiplet, which splits 
into 2Si @B 2 . 

According to the neutron diffraction experiments the 
magnetic moments lie predominantly in the ac plane 
with small component along b£2- Therefore, one can de- 
scribe the magnetic structure by two order parameters 
iVx,^,x) and (?7z,£z), which give the magnetic moments 
in the ac plane and transform according to B±. It can 
be argued that the observed small y magnetic moment 
component also results from the IR B±, which enters, 
however, into the other exchange multiplet Bi <S> V = 
Bi(B2B 2 - Both order parameters allow Lifshitz invariants 
rj/jd^/dx-^drj/s/dx and T)pd£p/dz-£pdrip/dz (/3 = x, 
z), which are responsible for the long- wavelength modu- 
lation in the ac plane. The magnetoelectric interaction 
is given by {r) x £,z — £xVz)Pp- The direction of the electric 
polarization appearing at T c was not determined exper- 
imentally yet, but it can be argued from the aforemen- 
tioned magnetoelectric interaction that the polarization 
should lie in the ac plane. 

Similar to NFSO the phase transitions in NFGO, which 
has the same crystal structure, can be described start- 
ing from the orthorhombic praphase C'mcm. In the 
k = (1/3,0,0) point of the Brillouin zone, which retains 
in the orthorhombic structure, the space group possesses 
four two-dimensional IR's Sj (i = 1,2,3,4). The per- 
mutational representation of the Fe 3+ ions is given by 
Si E4 and the magnetic structure is given by the ex- 
change multiplet £ 4 <g> V = £i®£ 2 ©£3- Thus, we define 
the order parameters (a x , b x ), (a y , by), and (a z , b z ) trans- 
forming according to IR's £3, £2, and Si, respectively. 
The fact that the order parameters (j) x , £ x ) and (rj Zl £ z ) in 
the monoclinic structure transform according to a single 



12 



IR is reflected in the existence of the invariant 

(a x a z + b x b z )U xz . 
The magnetoelectric interactions are given by 
(a x b z - b x a z )P z , 

(a x b z - b x a z )U xz P x . 



(10) 



(11) 



The description of the phase transition sequence in 
NFGO using the thermodynamic potential expansion re- 
sembles that of NFSO and will not be given here. How- 
ever, as follows from magnetoelectric invariants Eqs. (|10[) 
and (fTTj) the x component of the appearing electric po- 
larization P x is proportional to U xz , which may result 
in its smallness as compared to P z . The discussion on 
the influence of external magnetic field in the ac plane 
on the phase transitions and magnetoelectric properties 
given for NFSO is also valid in the case of NFGO. 



D. Cu 3 Nb 2 8 

The magnetoelectric effect was recently found in 
Cu3Nb2 0g, which possesses centrosymmetric triclinic 
structure described by the space group PI (Cj)£2- Upon 
lowering the temperature it undergoes two magnetic 
phase transitions at T/v ~ 26 K and T-i w 24 K lead- 
ing to two magnetically ordered phases, which we de- 
note AF2 and AF1, respectively. The neutron diffrac- 
tion studies reveal that the magnetic structure in the 
phase AF1 is incommensurate with the crystal lat- 
tice and can be described by the wave vector k m = 
(0.4876,0.2813,0.2029)^ The electric polarization ap- 
pears at Ti and reaches values about 17 /iC vaT 2 . The 
magnetic structure of the phase AF1 is characterized by 
Cu 2+ spins rotating in the plane approximately perpen- 
dicular to the reciprocal space direction (1,2,1)^ 

The phenomenological theory of phase transitions 
can be developed starting from the space group PI 
and using the closest commensurate wave vector k = 
(1/2,1/4,1/4). In this point of the Brillouin zone the 
space group possesses one two-dimensional IR GPi and 
the complete magnetic representation of the Cu 2+ ions 
located at positions Cui (la) and Cu2 (2i) is given by 
9GPi. In order to describe spins rotating in a plane one 
has to consider two order parameters (ai, b\) and (02, ^2), 
which transform according to the same IR GPi. The spa- 
tial modulation of magnetic order is due to the Lifshitz 
invariants aidbi/da — bidai/da (i = 1,2, a = x,y,z), 
whereas the magnetoelectric interaction is given by 

(ai& 2 - b 1 a 2 )P a - 

The triclinic structure can be transformed to the mon- 
oclinic one described by the space group C2/m (C^) 
as shown in the Table Q1J 34 The phase transition C2/m- 



TABLE II: The left panel shows the atomic positions of 
CU3ND2O8 and lattice parameters in the triclinic PI struc- 
ture, whereas the right panel gives those, which result in 
the monoclinic C2/m structure. The monoclinic struc- 
ture is characterized by the lattice parameters a— 9.26 A, 
b=6.8 A, c=5.18 A, /3=109.9°, and atomic positions Cui 
- (0,0,0), Cu 2 - (0.1545,0,0.5403), Nb - (0.5964,0,0.7783), 
Oi - (0.0511,0.2323,0.7036), O2 - (0.7868,0,0.7365), O3 - 
(0.3568,0,0.8265). 
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o„ 


0.3554 


0.7383 
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0.8188 




0.2833 


o 4 


0.1735 
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0.1735 
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Pi is described by the order parameter U transform- 
ing according to IR GM 2+ and representing deformation 
tensor components U xy and U yz . In the Cljm struc- 
ture the magnetic order can be described by the order 
parameters belonging to the (—1/2,0,1/2) point of the 
Brillouin zone, in which the space group possesses two 
two-dimensional IR's Bi and B2. The permutational rep- 
resentation of Cui and Cu2 atoms is given by Bi and 2Bi , 
respectively. The multiplets corresponding to IR Bi split 
into Bi (g) V r/ =Bi®2B2. The y components of spins trans- 
form according to Bi, whereas the x and z components 
according to B2 (here the x and y axes of the orthogonal 
set [x, y, z) are chosen along the a and b crystal axes, 
respectively). The neutron diffraction experiments^ al- 
lowed determining the direction of the normal to the spin 
rotation plane in the phase AF1. In the monoclinic struc- 
ture the normal makes angles of 23°, 73°, and 124° with 
the crystal axes a, b, and c, respectively. The precise de- 
termination of the would be normal (if the structure were 
monoclinic) is not possible due to the slight ambiguity in 
the mutual orientation of the triclinic and monoclinic lat- 
ices. However it can be argued that the normal lies close 
to the ac plane. Therefore, in the AF1 structure both 
IR's Bi and B2 condense, which describe the magnetic 
moment components along y and in the xz plane, respec- 
tively. 

Similar analysis of the electric polarization direction^ 
reveals the angles 13°, 82°, and 121° with the a, b, and 
c axes of the monoclinic cell, respectively. Therefore, all 
the polarization components P a , p,, and P c are nonzero 
with the largest being P a . 

We denote by (pi,Qi) and (^2,92) the magnetic order 
parameters transforming according to Bi and B2, respec- 
tively. The fact that (ai,&i) and (02,^2) transform ac- 
cording to the same IR is represented by the invariant 



(j>iq2 - qiP2)U, 



(12) 
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whereas the magneto-electric interactions are given by 

(P1P2 + qiq2)Pb, (13) 



(j>iP2 + qiq2)UP ai {piP2 + qiq2)UP c 



(14) 



Therefore the first phase transition at T/v is con- 
nected with the simultaneous condensation of (pi,<7i) 
and {p2,q2) so that {p\qi — q\P2) 7^ 0, which results in 
(j>iP2 + qiQz) = and the absence of electric polariza- 
tion. In polar coordinates (p\ = ricoscj), q± — rising, 
Pi = ri cos(^)+A0), q% = r2 sm(cf)+A(j))) this corresponds 
to the phase shift of A<fi = ±ir/2. The analysis of the 
magnetic representation shows that the resulting mag- 
netic structure is characterized by the sinusoidal modula- 
tion of Cui and Cu2 spins with all three spin components 
x, y, and z being nonzero. The second phase transition 
at T2 results in A</> 7^ ±7r/2 and the appearance of all of 
the polarization components due to Eqs. (jTHJ) and (Tl4"l) . 
The polarization direction is close to the crystal axis a 
as shown above and, therefore, the coefficient at the re- 
spective invariant in Eq. (|14|) is the largest. The analysis 
of the thermodynamic potential is similar to that per- 
formed for pyroxenes above (Sec. IIVCI) and will not be 
given here. 

The experimental study of CusM^Os in external mag- 
netic fields has not been performed yet. However, ac- 
cording to our model the products of magnetic field com- 
ponents H x H y and H y H z transform according to GM 2+ 
similar to U and allow the interactions {p\q2 — qiP2)H x H y 
and {piq2 — q\P2)H y H z , which therefore directly influence 
the coefficient at / = (pi<?2 — 91P2) in the thermodynamic 
potential expansion. The decrease of the absolute value 
of the coefficient at / favors shrinking of the phase AF2, 
whereas its increase results in the growth of the AF2 
temperature interval and shrinking of the phase AF1. 



E. a-CaCr 2 4 

a-CaCr204 was recently shown to become magneto- 
electric below the antiferromagnetic ordering temper- 
ature Tn = 43 K i 59 i 60 It possesses an orthorhombic 
crystal structure described by the space group Pmmn 
(i?2ft)) m which chromium ions form a distorted trian- 
gular latticed Below T/v an incommensurate magnetic 
order appears with the wave vector k — (0,0.3317,0) 
and magnetic moments rotating in the ac planei 59 ' The 
electric polarization appears discontinuously below Tjv, 
though the polycrystalline nature of the sample did not 
allow the determination of the polarization direction^ 

The neutron diffraction experiments show that the 
magnetic order is described by two distinct IR's Ai 
and A 2 of the space group^ 9 - The phenomenological 
model of phase transitions in a-CaCr204 was suggested 
in Refs. [59.60]. Indeed, the magnetic structure can be 
described by the order parameters (ai,6i) and (02,62) 



belonging to the wave vector k — (0, 1/3, 0) and trans- 
forming according to IR's Ai and A2, respectively. In 
the following we define the x, y, and z axes paral- 
lel to the a, b, and c crystal axes, respectively. The 
electric polarization along the y axis is given by P y oc 
(ai<22 + 61 62) (01^2 — 6102), i-e., proportional to the fourth 
order of the magnetic order parameters, which is the pos- 
sible reason of low values of polarization of the order of 
0.5 /iC/m 2 ^ 

Thus, two different IR's Ai and A2 condense at Tjv 
as a result of a weakly first order phase transition! 59 i 60 
However, such transition is only possible if the coeffi- 
cients at a\ + b\ and a 2 . + h\ in the thermodynamic poten- 
tial expansion pass through zero almost simultaneously, 
which is rarely observed for two distinct IR's. Indeed, one 
can show that this could be the case. ct-CaCr204 pos- 
sesses two different sets of chromium ions at positions 
(1/2, 1/2, 1/2) and (0.4932, 1/4, 0.0046). 61 The permuta- 
tional representation for both of them in the (0,1/3,0) 
point of the Brillouin zone is given by Ai © A2 © A3 © A4. 
For both sets of ions the observed magnetic structure is 
described by the A4 © V exchange multiplet, which ex- 
pands into A 4 ©y' = AiffiA 2 ffiA 3 . The IR's Ai, A 2 , and 
A3 describe the components of magnetic moments along 
the z, x, and y axes, respectively. Thus, Ai and A2 be- 
long to the same exchange multiplet, which explains the 
closeness of the above mentioned coefficients. 

Further analysis can be performed by noting that the 
orthorhombic structure of a-CaCr204 can be represented 
as a slightly distorted hexagonal lattice. The lattice pa- 
rameters and the atomic positions of the orthorhombic 
structure are given in Table HTT1 The structurally similar 
compound /?-SrRh2 04 crystallizes in the hexagonal struc- 
ture PQz/mmc (Dg h ) and has strontium ions disordered 
over two positions 2b and 2c. 62 Similar praphase can be 
introduced for a-CaCr204, which has a very weakly dis- 
torted hexagonal structure as shown in Table IIII1 Thus, 
small ion displacements towards positions of higher sym- 
metry and the introduction of Ca disorder similar to that 
of Sr in /3-SrRIi2C>4 lead to the PQ^/mmc praphase. In 
the hexagonal structure the six-fold axis is parallel to the 
orthorhombic a axis and we define the orthogonal X, Y, 
and Z axes parallel to the — z, y, and x axes respectively. 
It can be shown that the order-disorder phase transition 
P&s/mmc-Pmmn is described by the order parameter 
(111,1*2, U3, u&, 115, ug) [phase state (0,u,0,0,0,0)[ trans- 
forming according to the IR A4 of the P63 / mmc space 
group and belonging to the (1/4, 1/4, 0) point of the Bril- 
louin zone. 

The two inequivalent sets of four chromium ions in 
the orthorhombic cell stem from the position 2a of the 
hexagonal lattice, which is supported by equal magnetic 
moment magnitudes found for these two sets by neu- 
tron diffractioni^i In the hexagonal praphase the ap- 
pearing magnetic order is described by the star of wave 
vector (1/6,1/6,0) (A point of the Brillouin zone). In 
this point the permutational representation of chromium 
ions is given by Ai © A2, whereas the observed magnetic 
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TABLE III: Atomic positions of a-CaCr 2 4 in the orthorhombic Pmmn (a = 11.0579 A, 6 = 5.8239 A, and c = 5.0553 A) 
and hexagonal P6s/mmc [a = 11.0579 A, 6 = 5.8306 A, and c = 5.0495 A (orthorhombic setting)] structures. In the table 
the atomic positions in the hexagonal structure are given with respect to the orthorhombic lattice. In the hexagonal setting 
(a = 2.9153 A, c = 11.0579 A) the atoms are at Cai - (1/3,2/3,1/4) (occ.=0.25), Ca 2 - (0,0,1/4) (occ.=0.25), Cr - (0,0,0) 
(occ -1.00), and O - (2/3, 1/3,0.1) (occ.=1.00). 





Pmmn 


PQz/mmc 




X 


y 


z 


occ. 


X 


y 


z 


occ. 


Cai 


3/4 


3/4 


0.3512 


1.00 


3/4 
3/4 
1/4 


3/4 
1/4 
1/2 


1/3 
1/3 
1/6 


0.25 
0.25 
0.25 


Ca 2 


1/4 


3/4 


0.0385 


1.00 


1/4 
3/4 
1/4 


3/4 

3/4 








1/2 


0.25 
0.25 
0.25 


Cn 


1/2 


1/2 


1/2 


1.00 


1/2 


1/2 


1/2 


1.00 


Cr 2 


0.4932 


1/4 


0.0046 


1.00 


1/2 


1/4 





1.00 


Oi 


0.4022 


1/4 


0.3365 


1.00 


0.4 


1/4 


1/3 


1.00 


o 2 


0.5904 


1/4 


0.6825 


1.00 


0.6 


1/4 


2/3 


1.00 


o 3 


0.5989 


0.4996 


0.1665 


1.00 


0.6 


1/2 


1/6 


1.00 



structure is described by the A2 ® V multiplet, which 
expands into Ai ® A3 ® A4. The IR A3 describes the 
component of magnetic moments parallel to the six-fold 
axis, whereas Ai and A4 the components in the hexagonal 
plane XY. The observed magnetic structure is given by 
A 3 and A4. Therefore, the introduction of the praphase 
does not lead to merging of IR's and the magnetic or- 
der is still described by two distinct IR's. Nevertheless, 
in the praphase two sets of chromium ions present in 
the orthorhombic structure merge to give only two ions 
per unit hexagonal cell (compared two eight ions in the 
Pmmn structure). 

The hexagonal praphase allows establishing the hierar- 
chy of interactions in a-CaC^C^. Let us denote by rji, 6 
and Ci (i = 1 — 6) the magnetic order parameters trans- 
forming according to IR's A 3 , A4, and A 1; respectively. 
The components 2 and 5 of these order parameters be- 
long to the wave vectors, which correspond to (0, ±1/3, 0) 
in the orthorhombic lattice. The appearing wave vector 
(0, 1/3,0) (in the orthorhombic lattice) allows for inter- 
actions linear in magnetic field and of third order with 
respect to the magnetic order parameters such as, for 
example, (af — ia\b\)H x or (a 2 b-2 + 2a\b\a,2 — b 2 b 2 )H z . 
However, one can show that when considered from the 
praphase they become proportional to u 2 . The first of 
these interactions stems, for example, from the invariant 
Hz{{u\ - ul){H - 36$) + {u\ - ui)(g - 366 2 ) + (uj - 
M l)(£f — 36£f), where u 2 — u and u,i — for i ^ 2. 
These anisotropic interactions could lead to rich behav- 
ior in applied magnetic fields, but it can be argued that 
coefficients at these invariants are small since they are 
proportional to u 2 and the distortion of the hexagonal 
lattice does not lead, for example, to the splitting of the 
energies of states corresponding to different IR's (A3 and 
A4) in the exchange multiplet. 

The orthorhombic symmetry allows the invariant 



H x Hz{o-ib2 — o-2bi), which stems from the interaction 

H x H z (mU - r)i£i + 2 V2& ~ 27756 + %&> - %6) 

+V3H Y H z (m£i ~ ?7i6 + 7736 ~ ?7e6)- 

(15) 

This interaction is not proportional to u 2 and may ex- 
plain the experimentally observed suppression of elec- 
tric polarization in magnetic field in polycrystalline 
samples.— Moreover, the invariants following from the 
praphase reflect the hexagonal pseudosymmetry present 
in cv-CaCr204, which can be tested by experiments on 
single crystals. 

Thus, it can be stated that the hexagonal praphase 
symmetry of ev-CaC^CU reduces to orthorhombic due to 
the ordering of Ca 2+ ions. This ordering does not lead to 
the splitting of energies of states corresponding to A3 and 
A4 stemming from the same exchange multiplet. The sys- 
tem possesses a magnetic instability in the (1/6, 1/6, 0) 
point of the Brillouin zone (in the hexagonal lattice), 
which has six nonequivalent wave vectors. Due to the 
orthorhombic distortion the appearing wave structure is 
described by two wave vectors only, which correspond 
to the (0, 1/3,0) modulation vector in the orthorhombic 
lattice. However, one can expect other arms of the six- 
vector star to appear under applied magnetic fields due 
to the interaction Eq. (|15p . for example, or other invari- 
ants such as i?z(666 + 666 - 666 + 666)- The 
two other pairs correspond to the (0, 1/3, 1/2) Brillouin 
zone point of the orthorhombic lattice. 



F. FeTe 2 O s Br 

FeTe20sBr (FTOB) presents an interesting case of 
magnetoelectric material. 63 FTOB possesses a mono- 
clinic structure described by the space group P2\/c 
(^2/1 )• Upon lowering the temperature it experiences 
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two magnetic phase transitions at Tjvi=11 K and 
T/V2=10.5 K, which result in HT and LT incommensu- 
rate structures, respectively. Both of the magnetically 
ordered phases are described by the modulation vector 
(1/2, 0.463, 0)i 63 i 64 Ferroelectric polarization appears in 
the ac plane below Tni with the largest component along 
the c axis (F c =8.5 ^C/m 2 , P a =l ^C/m 2 )^ 

The modulation vector of the magnetically ordered 
phases is close to the commensurate value k — 
(1/2, 1/2,0), which we use for the description of phase 
transitions in FTOB. In this Brillouin zone point the 
space group P2\jc possesses a single two-dimensional IR 
C\ and we employ two magnetic order parameters (ai, 6i) 
and (a 2 , b 2 ) transforming according to C\. The order pa- 
rameters allow Lifshitz invariants 



I 



LI 



a i 



dbi da\ 



dy 



dy 



'L2 



a 2 



dy 



da 2 
' dy 



which are responsible for the long- wavelength modulation 
along the y axis. 

The magnetoelectric interactions are given by 

P y (aib 2 - ha 2 ) = -P y r\v 2 sin(0i - <f> 2 ), (16) 



P a {a\b 2 + ha 2 ) = P a nr 2 sin^x + 4> 2 ), (17) 

where a = X or z and we assumed a\ = r\ cos^i, b\ — 
T\ sin (j)\ , a 2 — r 2 cos <p 2 , and b 2 = r 2 sin <f> 2 . 

The expansion of the thermodynamic potential up to 
the fourth order in powers of the order parameters can 
be written in the form 



$ = 



1 f f at T a 2 . 

V / I 2 2 
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7lj2 , ft 
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(18) 



where ai, a 2 , k, /3i, /3 2 , 71, 72, /, 01, a 2 , Si, and 5 2 



are phenomenological coefficients, I\ = 



a\ 



bir 



bf, J 2 

^2 , 






h = 



a\a 2 + b\b 2 Ji = c 
J : (oifc + b x a 2 )\ l n = (9 ai /aj/) 2 + (96i/ai/) . 
J52 = (da 2 /dy) 2 + (db 2 /dy) 2 . The analysis of the ther- 
modynamic potential Eq. (|18[) shows that upon lowering 
the temperature the paramagnetic phase becomes unsta- 
ble with respect to simultaneous appearance of spatially 
modulated (ai,6i) and (a 2 ,b 2 ) with n 7^ 0, r 2 7^ 0, 
0i = k\y and c/>2 = k 2 y + A<j>, where k\ — k 2 . The ap- 
pearing magnetic structure is characterized by A</> = 
and can be associated with the HT phase. The magne- 
toelectric interactions Eqs. (TiTj]) and (fTTj) imply that this 
phase is paraelectric since Eq. (|16p is identically zero, 
whereas Eq. (1171) averages out to zero upon integration. 
Further temperature decrease may result in the phase 
transition to the phase with A</> 7^ and the appearance 
of P y . However, the experimental data reveals that the 
LT phase is ferroelectric with P x and P z nonzero. 65 The 



interaction Eq. (jTTj) implies that such electric polariza- 
tion is possible if k\ — ~k 2 and A</> 7^ 0. This magnet- 
ically ordered phase can be indeed realized if a\ and 172 
possess different signs. Taking into account that (ai,6i) 
and (a 2 ,b 2 ) belong to the same exchange multiplet we 
assume for simplicity ot\ = a 2 = a and o\ = — a 2 = a. 
Thus, the phase transition to the nonferroelectric HT 
phase n 7^ 0, r 2 7^ 0, and (j}\ — 4> 2 = ky occurs at 



with 



a c i 



k = 





CT 2 

25 
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+ 2a 2 


V(T 2 - 
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V25 



If / < 0, subsequent temperature decrease results in the 
phase transition to the ferroelectric LT phase with n 7^ 
0, r 2 7^ 0, 4>i = ky, and <fi 2 = —ky + tt/2, in which 
the electric polarization lies in the ac plane according to 
Eq. (|17l) as observed experimentally. In the limit of small 
k this phase transition occurs at 



a = a cA 



25A 2 k 2 ct 2 (A 1 +A 2 -16/) 



2A 2( 7 4 (A 2 - 8/) 2 - S 2 k, 2 (A 1 A 2 - 64/ 2 )(A 2 - 4/) ' 

where A ± = 3/3i + 47! and A 2 = 3/3 2 + 47 2 . The LT phase 
is characterized by k = —a/ (28). Therefore, the modula- 
tion vector experiences a jump at the HT-LT phase tran- 
sition of the order of |<5k 2 /(4<t 3 )| in the limit of small n. 
The HT-LT phase transition is, thus, a first order tran- 
sition, which is experimentally confirmed by the abrupt 
change of the magnetic structure at T^ 2 . 63 



V. DISCUSSION 

In Sec. HVI we have performed the representation anal- 
ysis of several magnetoelectrics and compared it with 
the available experimental data. In all the considered 
cases the electric polarization is induced by two mag- 
netic order parameters. In some cases (MnWC>4, CuO, 
and a-CaC^O-j) they transform according to different 
IR's of the crystal symmetry group of the paraphase, 
whereas in the remaining cases (pyroxenes, CusM^Os, 
and FeTe205Br) according to the common IR. The anal- 
ysis of the exchange Hamiltonian symmetry shows that 
in all the studied cases, as well as in CuCl2 2:L and magne- 
toelectrics delafossites, 6 - the relevant IR's describing the 
magnetic structures at zero applied external magnetic 
field belong to the same exchange multiplet. Therefore, 
according to the concept of a single irreducible represen- 
tation proposed by Landau^i in the considered magne- 
toelectrics the magnetic structures at zero field are de- 
scribed by a single IR though not of the space group of the 
paraphase, but of the symmetry group of the exchange 
Hamiltonian. This conclusion is in accordance with the 



16 



observation that such is the case in majority of magnetic 
crystals^ However, it should be noted that in roughly 
90% of magnetic crystals the magnetically ordered phases 
are described by a single IR of the space group of the 
crystal, which is a stronger statement and is a conse- 
quence of considerable splitting of exchange multiplets by 
anisotropic interactions^ Therefore, it can be concluded 
that it is an experimental fact that in magnetoelectrics 
the splitting of exchange multiplets by anisotropic inter- 
actions is small, which results in the closeness of insta- 
bilities with respect to different IR's entering into the 
exchange multiplets. 

The splitting of exchange multiplets is determined 
by anisotropic interactions, which possess the symme- 
try of the crystal^ Therefore, latent pseudosymmetries 
present in magnetoelectric crystals may have strong in- 
fluence on their magnetic properties. Earlier we have 
shown that the crystal structures of some magneto- 
electrics can be represented as slightly distorted struc- 
tures (praphases) of higher symmetr y 19 i 20 In this work 
we further developed this approach and used it in the de- 
scription of MnWC>4, CuO, pyroxenes, Cu3Nb 2 0s, and 
a-CaCr204. Introduction of the praphase has several ad- 
vantages. In the cases of M11WO4 and CuO two IR's 
of the space group of the paraphase, which describe 
the magnetic structures, merge into a single IR of the 
praphase space group. Therefore, the observed magnetic 
phase transitions in these magnetoelectrics can be inter- 
preted as between the phases induced by a single order 
parameter. In the cases of pyroxenes and Cu3Nb2 0s two 
relevant order parameters transform according to differ- 
ent IR's of the praphase space group, whereas they trans- 
form according to the same IR when considered from 
the paraphase. Thus, two different IR's of the praphase 
merge into a single IR upon a phase transition to the 
paraphase. Finally, in the case of a-CaCr2 04 two rele- 
vant IR's being different in the paraphase description re- 
main different when considered from the praphase. How- 
ever, in this case the description from the hexagonal 
praphase significantly reduces the number of chromium 
ions per unit cell. 

In all the cases except a-CaCr 2 04, the reduction of 
the crystal symmetry from the praphase to the observed 
paraphase is described by a tensor quantity U (a com- 
ponent of the deformation tensor), which directly influ- 



ences the interaction of the two magnetic order param- 
eters. Therefore, external electric or magnetic field or 
elastic stress, which possesses the same symmetry as U, 
introduces additional symmetry allowed contribution to 
U. The influence of such external effects on the phase 
transitions were studied in detail for MnW04 and CuO 
and magnetic fields in the ac plane with H X H Z ^ 0. 
It is found that the experimentally observed phase di- 
agram with the magnetic field parallel to the easy axis in 
MnW04 is well described qualitatively by the suggested 
model. The resulting magnetic field-induced phase tran- 
sitions are due to the change in the splitting of the IR P\ 
of the orthorhombic praphase. Similar mechanism is sug- 
gested for the magnetic field-induced phase transitions in 
CuO. This allowed building temperature-magnetic field 
phase diagram of CuO, which can be verified experimen- 
tally. 

In the case of a-CaCr204 the praphase-paraphase 
phase transition is described by a non-tensor order pa- 
rameter with k ^ 0. Nevertheless, the presence of 
the hexagonal pseudosymmetry allows drawing impor- 
tant conclusions concerning the behavior of a-CaCr204 
in external fields. 



VI. CONCLUSIONS 

In summary, we performed representation analysis of 
magnetic phase transitions in several magnetoelectrics 
using the praphase concept and taking into account the 
symmetry of the exchange Hamiltonian. This approach 
allowed explaining the closeness of successive magnetic 
instabilities on the thermodynamic path and clarifying 
the behavior of these magnetoelectrics in external mag- 
netic fields of certain directions. Phenomenological mod- 
els of phase transitions for some magnetoelectrics were 
suggested. 
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